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We study the finite temperature Casimir interaction between a cylinder and a plate using the exact 
formula derived from the Matsubara representation and the functional determinant representation. 
We consider the scalar field with Dirichlet and Neumann boundary conditions. The asymptotic 
expansions of the Casimir energy and the Casimir force when the separation a between the cylinder 
and the plate is small are derived. As in the zero temperature case, it is found that the leading 
terms of the Casimir energy and the Casimir force agree with those derived from the proximity 
force approximation when rT S> 1, where r is the radius of the cylinder. When aT <C 1 ^ rT (the 
' medium temperature region), the leading term of the Casimir energy is of order whereas for the 

Casimir force, it is of order Ta. In this case, the leading terms are independent of the separation 
' When 1 ^ aT <C rT (the high temperature region), the dominating terms of the Casimir energy 

^ \ and the Casimir force come from the zeroth Matsubara frequency. In this case, the leading terms 

are linear in T, but for the energy, it is inversely proportional to 02, whereas for the force, it is 
inversely proportional to a 2 . The first order corrections to the proximity force approximations in 
different temperature regions are computed using perturbation approach. In the zero temperature 
(-H ' case, the results agree with those derived in [Bordag, Phys. Rev. D 73, 125018 (2006)]. 
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I. INTRODUCTION 



In the past ten year, there has been a tremendous interest in the Casimir interaction between two objects. This is 
in part due to the advent in the precise measurement of the Casimir force [ll-Q. The breakthrough in the methods for 
' effective computations of the Casimir force has helped to make possible the exact and numerical computations of the 
, Casimir force beyond the precision afforded by the proximity force approximation. Some of the methods developed 
•■ for computing the Casimir interactions include the semi-classical approximation JX), 11], the optical path method 
'. [l^ - tl^ - the worldline approach [isl - lisl . the functional determinant or the multiple scattering method [l9l - l27i |. and 
. the exact mode summation method [28l [29| . These methods have been used to compute the Casimir interactions of 
a number of geo metric configurations. Among the most pop ularly studied configurations are the sphere-plane 
m, MM, MMiMi. the cyhnder-plane [S [ll, E [H H^M IH, ^ , the cylinder- cylinder [Hii M, 113, S| and the 



sphere- sphere 
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, ^ [2^ - ^61 itol . |40| configurations. It has been shown that all the different approaches give leading 
, y , order approximations to the Casimir force at zero temperature that agree with those obtained using proximity force 
rS ' approximations. There is an intense interest in computing the corrections to the proximity force approximations. 

[, Among the methods mentioned above, the functional determinant or multiple scattering approach and the mode 
" " " summation approach can give an exact formula for the Casimir energy or the Casimir force that is valid at all 
separations between the objects. In general, it is easy to deduce the asymptotic behavior of the Casimir force from 
the exact formula when the separation between the objects are large. When the separation is small, this become 
numerically complicated. Therefore, this hampers the comparison between the results obtained numerically using 
approximation methods and the result computed from the exact formula. In the pioneering work [To'l, Bordag has 
introduced an analytical method to obtain the first order correction to the proximity force approximation for the 
configuration of a cylinder in front of a plane. This is later generalized to the configuration of a sphere in front of a 
plane 

In recent years, there is a growing interest in the thermal correction to the Casimir interaction from both the 
theoretical and the experimental sides fl','?']. However, the finite temperature Casimir interaction between two compact 
objects have been less studied compared to the zero temperature case. Weber and Gies |41|:|4^] extended the worldline 
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algorithm to the finite temperature case and showed that there is an interesting interplay between temperature and 
geometry. Both the sphere-plane and the cylinder-plane configurations were considered. However, as in the zero 
ternperature case, so far the worldline method has only been developed for the case of Drichlet boundary conditions. 
In Esl Sij , the functional determinant representation of the finite temperature Casimir interaction between a sphere 
and a plane was derived which also took into account dielectric materials. The Casimir interactions at medium and 
large separations were obtained half analytically and half numerically. In [i^ , Bordag and Pirozhenko studied 
analytically the finite temperature Casimir interaction between a sphere and a plane. The asymptotic behaviors 
at small separation is computed to the leading orders and they are found to be coincide with the proximity force 
approximations at medium and high temperature j45j . 

Among the different two-object configurations, the sphere-plane configuration is the most studied since this has 
been extensively used in Casimir experiments [H, 0I- However, it has been proposed that the cylinder-plane setup 
has some of the advantages of both the parallel-plane and the sphere-plane setups [1, lU 13 • [S| , Klimchitskaya 
and Romero have explored the possibility of obtaining stronger constraints on Yukawa-type corrections to Newtonian 
gravity from measuring the Casimir force between a cylinder and a plane. Thus the cylinder-plane configuration is 
becoming increasingly important in the study of Casimir effect. The exact Casimir interaction between a cylinder 
and a plane at finite temperature has been discussed briefly in [21]. The low temperature asymptotic behavior of the 
thermal correction has been obtained for large separation. So far no work has considered the exact finite temperature 
Casimir interaction between a cylinder and a plane at small separation. The purpose of this paper is to fill in this 
gap. 

In this paper, we derive the functional determinant representation of the finite temperature Casimir interaction 
between a cylinder and a plane using the Matsubara formalism. The exact expression allows us to conclude that 
the Casimir force is attractive at any temperature for both the Dirichlet and the Neumann boundary conditions, 
and hence for the perfect electric conductor condition. We then derive the small separation asymptotic behaviors of 
the Casimir force in different temperature regions. The leading order terms are compared to those obtained using 
proximity force approximation. The first order corrections are also computed analytically. 

In this paper, we use units with ?i = c = feg = 1. 



In this section, we use the proximity force approximation (PFA) to find the leading asymptotic behavior of the 
Casimir force between a cylinder and a plate when the separation between them is small. This has been considered 
in [l^l but we use a different approach here that would show explicitly the leading contributions in each of the 
temperature regions. 



II. PROXIMITY FORCE APPROXIMATION 




J 



FIG. 1: The cross section of a cylinder in front of a plate. 



As shown in Fig. [1] a cylinder length L and radius r is placed in front of a plate. The distance between the cylinder 
and the plate is a, and the distance from the center of the cylinder to the plate is equal to a -|- r, which we denote by 
H. 
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The finite temperature Casimir force density for a pair of parallel Dirichlet /Neumann plates separated by a distance 
d is given by |50l] : 



480d4 



90 



k=l 1 = 1 



T:kl 

It 



(1) 



fc=l /=1 ^ 



(2) 



Eq. (HI is the low temperature expansion. It shows that in the low temperature region (i.e., dT <C 1), the Casimir 
force density is dominated by the zero temperature term 



The leading term of the thermal correction is the d-independent term 



90 



(3) 



(4) 



and the remaining terms go to zero exponentially fast as dT — )■ 0. Eq. ([3]) is the high temperature expansion. It 
shows that in the high temperature region (i.e., dT ^ 1), the Casimir force is dominated by the term 



which is called the classical term. The remaining terms go to zero exponentially fast when dT oo. 



(5) 




FIG. 2: The cross section of the cylinder-plate system in the xz-plane. 



Applying the proximity force approximation to the cylinder-plate configuration with either Dirichlet or Neumann 
boundary conditions, we need to integrate the Casimir force density between the parallel plates J^c^{a) over the 
rectangle {{x,y) : —r < x < r, < y < L}. At the point {x,y), the distance between the cylinder and the plate is 
a + r(l — COS0), where x = rsinO (see Fig. [5]). Therefore, the proximity force approximation for the force is 



Let i = 1 - cos 61. Then 



pPFA 



2L 



•^Cas(a + ''(l-cos0))(ix. 



dt 



V2t - t2 



(6) 
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For the Casimir energy, the proximity force approximation gives 

dt 



V2t - ' 

where Sq^^ is the finite temperature Casimir energy density for a pair of paraUel Dirichlet/Neumann plates given by 

1=1 k=l ^ ^ 



oo oo 



Obviously, 



oil - ^ T f nlT 1 



QTTiPFA 

pPFA _ "^Ca.s 
^Cas — 



e 



da 

As in [4^, we study the proximity force approximation in the following three regions: 

1. Low temperature: aT <C rT <C 1, 

2. Medium temperature: aT <C 1 <C rT, 

3. High temperature: 1 <^ aT <^ rT. 

A. The low temperature region 

In the low temperature region, the Casimir force is dominated by the zero temperature term. Substitute the zero 
temperature Casimir force density ^ into (jH]), we find that the proximity force approximation of the zero temperature 
Casimir force is 

pT=o,PFA ^ _ TT^Lr 1 _ dt 

240 Jo (a + rt)4^ ' V2t^' 

Making a change of variables 1 1-^ au, we find that 

^T=o,PFA ^ _ TT^Xm 1 1-au 



240 Jo {a + aru)4 ^2au ~ a^u^ 
' zdu — 



240\/2a2 Jo (1 + ru)40I TGSV^ai ' 

This agrees with the leading term found by using exact method in [19]. For the thermal correction, the leading term 
in the low temperature region comes from the distance independent term ([¥]), which gives 

For the remaining terms in ([!]), since aT <C rT <C 1, 

»»(-(^)*' 

for any < i < 1. Therefore, they give exponentially small contributions. In other words, we find that in the low 
temperature region, the proximity force approximation gives 



768%/2a5 



fEI^ = -^:::^ + .... (8) 
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and the leading term of the thermal correction is 

AtFEI^^-^LZIL + ,.,_ (9) 
Notice that this term is independent of a. In the same way, we find that for the Casimir energy, 

^''^ 1920V2ai 
For the thermal correction to the Casimir energy, we have 

The first term is independent of a and so it does not contribute to the Casimir force. The second term gives rise to 
the leading contribution to the thermal correction of the Casimir force 



B. The medium temperature region 

In the medium temperature region, one cannot ignore the contribution to the proximity force from the exponential 
terms in ([1]). In fact, when aT ^ 1 ^ rT, the function 



exp 



(a + rt)T 



is small if t ^ but it becomes ~ 1 when i — > 1. The contribution of these exponential terms to the proximity force 
approximation is given by 

^ ^ fc2 1 / nkl \ ^ , dt 
/=7rLrr> > — / -^exp — 1 - t , 

a^oirLT^^k^ 1 ( Txkl\ , , dt 



EEt 1 73-p h;:?^ (1-^) 



^ Jo "V rTtJ^ 'V2t^ 



ttLT ^ ^ fc2 / ttM \ t{t - 1) , 



fc=i 1=1 

To find the asymptotic behavior of this term when rT ^ 1, we use the inverse Mellin transform formula 
which gives 



C+'tOO 

r(z)u-^dz, (10) 



where A{z) is the function 



1=-^— r{z)( — )(:n{z-2)CR{z + l)A{z)dz, (11) 

J C— ZOO \ ^ / 



Aiz)^l^ -J^dt. 



A change of variables iH'(t + l)/2 gives 

/■°° dt 

A{z)^r-'j^ it + iy-^-{t-i)— 

^2^-3 ( r dt - r-jzi—dt+ [\t+\Y~^{i-t)% 

[Jo it + ir-' Jo it + iy-' Jo^ ' ^ wt, 



^2-(-V^(.-3)L|^+/V + l)-(l-.)^^ 
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The last integral is convergent for any z. Therefore, Aiz) has poles at ^ = 5/2,3/2,1/2,-1/2,.... Moreover, for 
J = 0,1,2,..., 

H«„,_,(r,=MW,.^H)i(i±i)(i-). 

The integrand in PT|) has simple poles at z = 5/2,3/2, 1/2,-1/2,... and at z = and z — 3. Evaluate the residues 
gives the asymptotic behavior for / when rT ^ 1. Namely: 



/ M C»(4).4(3) + C« (-2) Ito r(=)yl(j) 



j=0 ^ ' 

Since A{3) = 1 and Ci?.(-2) = 0, we find that 



2^ + 5 J 



The first term, — — — , cancel with the term Q coming from (jlj. Therefore we find that if aT <C 1 ^ rT, the 
asymptotic behavior of the thermal correction is 

..^v.^c. (i) . (I) - ^^c. (4) (I) - (4) c. (H 

In particular, the leading term is 

AtFZ^ ~ LTiV^^Cn (0 ~ -0.2778LTiV^. (12) 

This agree with the result obtained in [l^l . The leading term of the Casimir force is still given by the zero temperature 
term @ since aT < 1. 

For the Casimir energy, a similar computation computation gives 

,PFA Ly^Ti /5V /3\ iLr-^Ti 



which is independent of a. In particular, the leading term of the thermal correction to the Casimir energy is 

pPFA _ L^Ti ^ f5\ ^ f 3 



^tE^Is - -^f^C/?, j Ci^ ) = -0.1972LT*V^. (14) 



C. High temperature region 

In the high temperature region, the Casimir force is dominated by the classical term. Substitute the term ([S]) into 
)]), the proximity force approximation gives 

PPA Cfl(3)LrT 1 dt 



Making a change of variable f 1— >■ au, we obtain 



47r 7o (a + rau)^ y/2au - a^v? 

CR{i)LrT f°° 1 3C/^(3)i^/F^ 



4^/27^ai Jo {l + ruf^/^ 32^/2ai 
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For the remaining terms in ([2]), since 

exp (-47r(a + rt)T) 

is exponentially small when 1 <C aT <^ rT, they give exponentially small contribution to the proximity force approx- 
imation. In other words, in the high temperature region, the proximity force approximation gives 



For the Casimir energy, one has 

z?PFA 

The results obtained above are for a cylinder above a plate with Dirichlet or Neumann boundary conditions on 
both the cylinder and the plate. For perfectly conducting cylinder and plate, all the results have to be multiplied by 
two due to the TE and TM modes. 

III. EXACT FINITE TEMPERATURE CASIMIR ENERGY 

In this section, we consider exact expressions for the Casimir energy of the cylinder-plate system. The exact zero 
temperature Casimir interaction energy of the cylinder-plate system is given by 

^caf =2^j^ Tr In (l - M ( VC^+P) ) d^dk ^ — ^Tr In (1 - M{0) d^ (16) 

where M(^) is an infinite matrix with elements Mjk{^), — oo < j,k < oo. 



for Dirichlet boundary conditions, and 



for Neumann boundary conditions. For perfectly conducting boundary conditions, the Casimir energy is the sum 
of the Casimir energy for Dirichlet boundary conditions and the Casimir energy for Neumann boundary conditions. 
These formulas have been proved by using the mode summation approach 29] and the path integral approach [19, 21j] . 
Using the Matsubara formalism, one then finds that the finite temperature Casimir energy is given by 



V £' 1^ (l - ^ ( + fc") ) dk 

° (17) 
TT 1'°° C 



Ecas — 

TT 

1=0 



where = 2t:IT are the Matsubara frequencies. This formula is suitable for studying the high temperature behavior 
of the Casimir energy. In particular, in the high temperature region, the Casimir energy is dominated by the classical 
term corresponding to zero Matsubara frequency: 

Ec.s -ES^s = ^ Tr ln(l - M(0)de 
To study the low temperature behavior, one can use the Poisson summation formula 

oo |.oo ^ poo 

= / f{x)dx + 2^ f{x)cos2Trlxdx 
1=0 1=1 
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with 



= / ^7r4^=^Tr ln(l - M{0)d^. 
J27VXT - (2tixTY 



/>oo 

Then the term / f{x)dx reproduces the zero temperature CasnTiir energy ()16p . The thermal correction to the 
Jo 

Casimir energy is 

AtEc.s = 2 f{x)coa2TTlxdx^-^ / ^=L=Tr ln(l - M(OKcos -da: 
Jo 2t:^T J„ _ 3,2 T 

I e I ^_ cos |d.Tr h.(l - MiO)d^ (18) 

Therefore, the low temperature expansion for the Casimir energy is given by 

Ec.s = + ^ E ^-^0 (I) Tr ln(l - M(0)de 

Alternatively, one can apply the Abel-Plana summation formula to the first line of (|17l) to obtain 



^ .00 /.oo i (Tr ln(l - M(iv/|^^^)) -Tr ln(l -M(-iV^^^^)) 



i?Cas =i?Jaf + 7^ I I ^ ^ ^ ^d^dk. 

Jo Jk er — 1 

The second term is the thermal correction Ayi^cas- By a change of variables, it is equal to 
Since 



we find that 

L 



^^"{*[Tr ln(l - M(zO) - Tr ln(l - M(-zO)] }^A'o ^C- 



Notice that this term can be obtained from ([TS]) by a rotation ^ — > «^ using Jo(2) = Rei/p (2) and [iz) = 
2i 

Ko{z). The integrand does not have singularities since we are considering an open geometry whose spectrum is 

purely continuous [45|. 

For the Casimir force, one has 

= - ^ £' r f ^ - ««))-'Q«)} «. 



where 



for Dirichlet boundary conditions, and 

QMO - QUO = §^^K',+k{2m) 



for Neumann boundary conditions. Since Mjk{^) and QjkiO are positive functions, using (1— M) ^ = 1+M+M^ 
one concludes that the Casimir force is always attractive at any temperature. 
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IV. ASYMPTOTIC BEHAVIOR OF THE CASIMIR ENERGY AND CASIMIR FORCE AT SMALL 

SEPARATION 



Define the dimensionless parameter 



a 
r 



The first two leading terms of the zero temperature Casimir energy at small separation £ <JC 1 has been obtained by 
Bordag in [l^. It was found that for Dirichlet conditions on both the cylinder and the plate, 



E. 



D,T=0 



36 



whereas for Neumann conditions on both the cylinder and the plate 

,Ar,T=o 



E. 



Cas 



For the Casimir force, one then obtains 



F, 



F, 



D,T=Q 
Cas 

N,T=0 
Cas 



1920V2ai 

768^/2ai 
768%/2at 



1 



7 40 
36 ~ 3^ 



(19) 



(20) 



60 
T_ 
60 



In the following, we would derive the asymptotic behavior of the Casimir energy and the Casimir force at finite 
temperature. This has been studied numerically in [42| using worldline method. The approach used here is similar to 
the one used in [l^ for the zero temperature case. The leading asymptotic behavior would be computed analytically 
using the exact formula for the Casimir energy and compared to the proximity force approximation. We would also 
discuss the limitation of this method when trying to obtain the first order correction analytically. 

The starting point is the formula (|17l) . Making a change of variable ^ = ui/r, one obtains 



-C/Cas — > ^ 

■rrr ' ^ 



1=0 



:Tr ln(l - A(cj))dw, 



where 



WT^^-+fc(2c.(l + e)), 



K'Aru) 



K,+k{2Lu{l + e)). 



Expanding the logarithm, we find that 

T,r ^S + l^l^. ./„,2, 



72^ ■ ^ 



Debye asymptotic expansions of Bessel functions give an expansion of Ai,-^iy.-^{uj) in the form: 



J [uj)duj. 



1 



r) exp(- 



271^ V2Vz \ (1 + a;2)2(l + cjs)^ 

Ui{t{uJi)) Ui(t(w2)) Ui{t{uJz)) 



X 1 + 



V2 



^3 



(21) 



V2 



Vl[t{LOl)) Vl{t{uj2)) Ul(t(w3)) 



cxp(^-i^3r;(w3) + J^i?7(a;i) + ^211(^2)^ 



X 1 + 



1^1 



1^2 



V3 



(22) 
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where 



bj u! 2(1 + e)uj 

J/1 1^2 1^3 



z , . 1 



7(z) ^^/l + z^ + log 7==, t{z) 



l + Vl + Z^' Vl + z2' 



"iW = — ^1 — ' 



24 ' ' ' 24 

Observe that after interchanging the subscripts 1 and 2 in the first fine of (|22l) , it coincides with the first fine of (|2ip . 

One can show that the expression —v^rjltu^) + virj^ioi) + i'2i]{'-^2) has a maximum value of when vi ~ V2 and £ = 0. 
This suggest that we rename jo as m, and as m + for 1 < i < s, and treat ni, . . . , as perturbed variables. 
Making a change of variables w i— >■ mw, and replacing the summation over nt to integration (which is the leading term 
after applying Poisson resummation) , we find that 

^Cas —— Y,' Y '"^ L, / / ■ ■ • / n ^m+n.,m+«.+ i (maj)d71i . . . dUsdoJ. (23) 

s=0 /=0 m=0 "^-S- i/a;2 ^ J -00 -00 

V 

Treating e and tii, . . . ,71^ as perturbed variables, one can expand A„i+„. „i+„.^j (mcj) with the help of computer. In 
fact, for the term — 1^377(^3) + t^i 77(^1) + i'2'n{^2), where ui = to + n^, 1^2 = to + nj+i, 1/3 = 2m + + we have an 
expansion of the form 

Ni N2 
k=0 j=0 

where ^kji^i, rii+i) is a homogeneous polynomial of degree k in rii and rii+i. The leading terms are 

4tov 1 + 1^ 

From here, we see that the leading contribution comes from m ^ e~^, ni,ni+i ~ m^/^ ~ e"-'^/^. Consequently, the 
term £^m^~^^kj{ni, n^+i) is of order £^^^~^ . To keep everything to the order e, we need to expand up to the terms 
with j + k/2 < 2. Doing similar order analysis to the other terms, and expanding the terms that are of order y/e and 
e in the exponential using — 1 + a: + a;^/2 + . . ., we find that 

Ajn+ni,ni+ni+i{inLu) ^ — —— ^ exp -2eTO\/l + - ' ' =r (1 + c{ni,ni+i) + b{ni,ni+i) + ...), 

2\/7rm(l + w^j 4 \ 4TOVl+t^ / 

(24) 

where c{ni, rii+i) and b{ni, ni+i) are terms of order ^/e and e respectively. In the case of Dirichlet boundary conditions, 

Dr ^ rii-Sn^+i {rii + ni+i){ni ~ rii+i)'^ 1 

c n,,ni+i)=— — — — + —-77; ;-3 £ , ^ (rir+ni+i) 

4TO(l+a;^) 8to2(1+w2)2 VI + 



6^(ni,ni+i) 



32m2(l+cj2)2 2(l + a;2) 

(nj - 3ni+i)(ni + nj+i)(ni - Ui+iY (n^ - 3ni+l)(7^i + rii+i) 

- 5 £ 3 

32m3(l+cj2)2 4to(1 + w2)2 

" 128TO4(l+a;2)3 +^ 2(TTi:;^("^ + "^^^^ ' " 8to2(1 + 0.2)2 + (r.. - n.+O 

2(71^+^1+1)^ 2\ ^ + 10rii7^.i+i + 7nf+i) 2 "1 2 - 3i:j^ 



4m(l+a;2)2 192to3(1 + a;2)f ^/\+UJ^ 16to(1 + 0.2)2 

Only the last term of b^{ni, Ui^i) comes from the second line of (|2ip . Therefore, for the case of Neumann boundary 
conditions, a comparison between ((2T|) and (|22l) shows that 



c^{ni,ni+i) = c^{ni+i,ni), b'^ (ni^m+i) ^ {m+i^rii) —j. (25) 

7n(l + o;^) 2 
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These relations between and c^, and between and can help us to derive the asymptotic behavior for the 
Neumann case from the asymptotic behavior for the Dirichlet case. 
Substituting ^ into we find that 



— Et ^E'E'™^ /, / (l+a.r^cxp(-2e(,s + l)myrT^)(l + ^-)d^, (26) 



Cas 



where 



^^-E^^^+E E 



s— 1 s 

IJ 1 



2«7r2TO2 (1 + w^)* 7-00 J-oo \ 4m\/r+7j2 



\2 



exp . b{ni,ni+i)dni . . . dus 



=lzr^-^77-—^ / ■ • • / exp c(?ii,n,+i)c(nj,nj+i)dni . . . dn. 



We have used the fact that 



— OO J —OO 



exp I ^"'-"-+^^''1 dn,... dn. ^'''''^^^ + "'^^ 



4my/l + ijp- ) \Js + 1 



.e/** is the term of order e. The term with order ^fe has dropped out since c(ni, n^+i) is odd in rii or n^+i. 

^1 and "^jj- are Gaussian integrals. They can be computed in the same way as explained in [l^. With the help of 
a machine, we find that for the case of Dirichlet boundary conditions, 

D e( ^, (2i + l)2\ 1 / 8i + 4 2422 + 242 + 7\ 1 2 ( ^ (2« + l)^\ 1 

' - ^-(4» + 2)+ ^ _ ^M +— -— — + +e^m 4^ + 2-^ ^ 



2V s + 1 y 16m V .5 + 1 (s + l)2 / (l + a;2)^ V «+l 7(1+^2)^ 

e/. -4i2 + 4i + 2 3(2i + l)2\ 1 1 /8i + 4 -24^^ + 12^ + 6 \h{2i + \f\ 1 



2V s + 1 (s + l)2yi + a;2 I6771 V s + 1 (s + 1)^ (s + 1)3 / (1 + a;2) 

sD 2 /^o- . 8ij+4i + 47 + 2\ 1 /Si + 4 12ij + 6i + 6j + 3\ 1 

' =e m 8i + 4 ' 



s + 1 /(l+u;2)i + 1 (s + l)2 y l + w2 

1 /36i + 18 60ij + 30i + 30j + 15\ 1 



^8m\{s + lf (.s + l)3 y'(l+a;2) 

Therefore, 

-2 

(l + a;2)^ 6V ' Il + oj^ 16m V s + iy(l + a;2)| 



f (. + If + 2{s + 1)) + Us + 1)^ + 2) ^ + ^ (-7[s + + 

V ' ^ V(i + a;2)i eV^ ^ yi + a;2 16m V 5 + 1/(1 + 0; 



(l+t^2) 



2H 



For the case of Neumann boundary conditions, the relations (P5|) suggest that we can make a change of variables 



rii ^ Tis+i-i. Then it is easy to see that '^l-^ — ^^s^j s-i ^'^d 



■,N ^s.D Vs + l [°° f°° ( {ni-n,+if\ ~ J J 

^ 'exp — , — T-rdni ...dus 



s—% 



2''7r2m2(l +o;2)4 \ 4mVl + ^2 / m(l + u;2)f m(l+w2)! 

Therefore, 

^sM _ ^s,D ^ ^ _(^ ^ 

m(l + (u;2) 2 

Return to the Casimir energy (|26p. the inverse Mellin transform formula gives 

/ \ 1 /-c+'oo 

exp(-2em(s + l)Vl + w2 =7^ / r(z)(2e)"^m"^(s + 1)^^(1 + w^)^ 2 rf;^. 

V / 2711 Jc-ioo 
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Therefore, 
exp (^-2£mVTTw2) A''° 

= - / 



c+ioo 

r{z + 2)(2e)-^-2m-"-2(s + l)-^-2(l + uj^^ 

C—lOO 

c-\-ioo 



((, + 1)3 + 2(. + 1)) — 1-^ 

^ (1 + ^ )^ 



dz 



c—too 
c+icx) 



dz 



r(2)(2£)-^m-^(s + 1)-^(1 + w^)-* 
1 



16rn 



7(s+l) + 



1 1 , . X 



+ 1/ (l+w2)i 16m 



s + lj (l+a;2)l 



1 

'2^ 



C+iOO 



r(0)(2£)-^m-^-i.^ (l+w^)-^ 



(2^ + 7)(2^-3) 1 _^ (2£ + l)(20 + 3) 1 



+ (l+a;^)-^ 



48 (s + 1)^-1 

" 2z-3 1 



24 (s + 1)^+1 

4z + 3 1 



12 (s + 1)^-1 48 (s + 1) 



>dz. 



Using 



we finally obtain 



1 . 0F m^M-i r(/i-i) 



f 

V TO^ 



where -Ecag is the leading order term and E^^^ is the first order correction term. For the case of Dirichlet boundary 
conditions, they are given respectively by 



TL^ 1 



rpO,D _ ^ ^ \ ^ 

Cas ~ 27rr ^ 2-Ki 



r(.)(2£)- ^^^"^ 



■dz, 



and 



E. 



1,D 
Cas 



rp T ^ -t />c+ioo 
„_n ' Jc-loo 



2-Kr 



s=0 



r(^) 



r(^) (s + 1)^+1 



(2z + 7j(2--3) 1 _ (2z+l)i2r + 3) 1 

1 "r 



(27) 



24 



+ 



48 (s + 1) 

+ r_2z-3 1 
r(^) [ 12 (s + 1)^+5 ~ 48 {s + iy+i 



{s + l) 



z+i 



Az + 'i 1 



>dz, 



(28) 



where 



&i{z)=V 



2z- 1 



00 00 

E'E' 



/=0 m= 



(m2 + r2^f) 



(/,m)/(0,0) 
00 00 



^2 



i=0 m=l (™2 +r2ef) 4 



For the case of Neumann boundary conditions, the leading order term is the same as the Dirichlet case, i.e., -Ecas ~ 



^Cas ' whereas for the first order correction term. 



^l,N _ ^l,D _ J_ 
Cas Cas 2-Kr ^ 

s=0 



c+ioo 



V{z){2e)- 



Mz) &i{z + l) 



(5 + 1) 



-dz. 



(29) 
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The functions ^i(z) and &2{z) have expansions (Chowla-Selberg formulas) of the form 



®l:l,,o(z) @'i:t.i(2)('"=0 term) 

2.-3 , , (30) 



S*1;I,.2(Z) 



2g + 3 



1 /2z + 3\ / 1\ 1 1 / I \~ flm 



^i;T,o(2)(/-0 term) ^1;T.i(2) (32) 

00 00 2s-3 



i=l m=l 



. f?i±l) c« (. + iV#r (?i^) c« (« - i I (2..r,.- 



®2,T.o(z)(i=0 term) ®2;r,l(2) (33) 

00 00 2z-l 00 00 2z-l 

v^E E (^) ^ (4-^'-?^) - E E ^^i;^^-^ (4-^^-^^) ■ 



l—l m—l l—l m- 



From these, one can deduce that &i{z) only has simple poles at z = 5/2 and z = 1/2, whereas S^2{z) only has simple 
pole at z = 3/2. 

Now we can derive the asymptotic behaviors of the Casimir energy at small separation in different temperature 
regions. First we consider the leading order term. Substituting ([50]) into (l?fl) . we can obtain the low temperature 
leading behavior. In fact, the first term in pop . ^^i-Lfliz), is the term corresponding to the zero temperature Casimir 
energy. The residue at the largest pole at z = 5/2 gives the zero temperature leading term: 



2^^^^ ^-2 \^ ^ ^ r(^) (s + 1)^+1 y 1920V2r2ef' 

which is exactly the leading term obtained by using proximity force approximation. Although i^i-L,a{z) also has a 
pole at z = 3/2, but this is canceled out by the pole at z = 3/2 of the second term i^i;L.i(z), which is the term 
corresponding to m = 0. 

Next, consider the contribution of the term &i-L,2iz) in ([501 to the thermal correction of the Casimir energy. Since 
this term does not have any poles, we find that the leading term of the thermal correction is of order e°, coming from 
the pole of r(z) at z = 0. This term would not contribute to the Casimir force since it is independent of e. Therefore, 
let us also consider the term of higher order in e by taking the residue atz = — 1,— 2,.... We have 



AT..i?8..--^EEi^-^--. r(z)(2s^-^^^^-(^) ' 



2-,t^^ '\ ' r(^) (. + 1)^+1 



2z + l^ 

2j + 3 



^ (-l)^2^-VLCfl(|-j) ^^ / I (lm\ (34) 



order e-^ term 
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Notice that for the residue at z = —j, j < 1, the summation over s is divergent, but we analyticaUy continued it to 
obtain the value ^^{(3/2 — j). 

The terms in p4|l are exponentiaUy small when rT <C 1. But we cannot use this to conclude that the thermal 
correction is exponentially small when rT <^ 1. The Debye asymptotic expansion cannot be used to study the asymp- 
totic behavior of the thermal correction in the low temperature region since all the terms in the asymptotic expansion 
would contribute, and summing over infinitely exponentially small terms would lead to polynomial order terms. Nev- 
ertheless, we can use p4p to study the leading behavior of the thermal correction in the medium temperature region 
where aT <C 1 ^ rT. Using the formula pUj) again, we have 

°° °° / I \ — ^^4^ / lm\ 



, 2mrT J * \ rT J 
1=1 m=l ^ ^ ^ ^ 

i=lm=l"''J \ ^ ^ 

OO OO pC+ioO /r) rp\2w 



tm^ dt 



(35) 



^1 m^l 



Notice that from the term we have 



2-,t^^ ^- ' r(i^)(. + i)^+i; 



where we have used the identity 



r(£)c.W = .--r(i^ 



This term cancel with the second term in the last line of psp . Therefore, we find that in the medium temperature 
region, the leading asymptotic behavior of the thermal correction is given by 



AV27: \2j^ \2j aV2tt \2 ) \2 J 8V2n \ 2 J \2 



(37) 



coming from the first term of ([55)) . In fact, one can also obtain these terms by using the term ^i;t,i in p2p and 
taking the residue of (P71) at j = 0, —1, —2, . . .. The first term in ([57]) is the leading term of the thermal correction 
to the Casimir energy in the medium temperature region. It agrees with the result of proximity force approximation 
(ITil) . For the Casimir force, taking derivative of ([57)) with respect to a gives 

Again, the leading (first) term agree with the result of proximity force approximation (I12p . In fact, the second term 
which is equal to O.lSSler^ri can also be predicted by proximity force approximation (see p^). 
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Now let us consider the high temperature region. In this case, the Casimir energy is dominated by the terms with 
I — 0. Substitute the / = term &i-t,o in (1321) into (l?fl) . and taking the residue at the largest pole of ^i;t,o at 
z — 3/2, we find that the leading term of the Casimir energy in the high temperature region is given by 



E 



0,cl 
Cas 



TL 
2'Kr 



5^Res,^3 r(z)(2e) 



l;T,o(^) 



1 



LT 



r(2^) (s + ir+i 



16\/2rei 



Ci?(3), 



agreeing with the result of proximity force approximation (fT5|) . As in the case of f^i;L,o, the pole of ^i;t,o at z = 1/2 
would actually cancel with the m = term which is not explicit in this case. 

From these discussions, we find that the leading order terms of the Casimir energy (P7)) which is derived from the 
exact formula agree with the proximity force approximations in the zero temperature, medium temperature and high 
temperature regions. Let us now consider the first order approximation from the term (1281) . For the zero temperature 
Casimir energy, substitute &i-^l,o{z + 1) and ^2;L,o(^) into (pS)) . and taking the residue at z = 3/2 give 



E. 



1,D,T=0 
Cas 



TL 



^Res,^Jr(z)(2e)- 



s=0 



h,L,o{z + l) 



(2z + 7)(2z-3) 



1 



(2z + l)(2z + 3) 



48 
2z- 3 



(s + 1)^+5 

1 4z + 3 



24 



12 [s + iy+i 



48 (s + 1) 



7 

X — e 



69120%/2r2£i imO^/ir^ei 36 



(39) 



for the case of Dirichlet boundary conditions. For the case of Neumann boundary conditions, ()29|) gives 



j^l,N,T=0 pl,D,T=0 



TL 
27rr 



^Res,= Jr(z)(2e) 



h-L.oiz) ^l;L,o(z + l) 



1 



s=0 



nL 



69120%/2; 



144 V2; 



1920^2; 



r(2^) r(2^) J(s + i)^+^ 

■K^L ( 1 40 

36 ^ 3^ 



(40) 



([39| and (l4pl) agree with the first order correction for the zero temperature Casimir energies (IT9l) and (l20t obtained 
by Bordag |19| . 

For the thermal correction, we have 



TL 



OO CX3 



j=0 s=0 



^2;T,l(^) 

r(2^) 



(2z + 7)(2z-3) 



1 



^i;T,i(2; + 1) 



24\/2i 



— C-R 



2z- 3 



48 



1 



12 



6\/2 



(s + 1)^+5 

4z + 3 1 
48 (s + 1)^+1 



(2z + l)(2z + 3) 1 

24 (s + 1)^+1 



for the case of Dirichlet boundary conditions. For Neumann boundary conditions, (j29p gives 



CX3 OO 



<'2;T,l 



(Z) ^l;T,l(2 + l) 



24^2) 



j=0 s=0 

1 

2 



18Cfl, 



r(^) r(^^ 

TxLrWi ^ ( 3 



(^ + 1) 



These terms are of order l/(rT) smaller than the leading term. For the Casimir force, we find that 



AtF, 



1,D 
Cas 



TxLTi 



= -0.1975Lr"^r5 
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and 



= -1.4290Lr-5r2 



For the high temperature (classical) term, substitute f^i;T,o(^ + 1) and i^2:T.o(^) into and taking the residue 
at ^; = 1/2 give 



E, 



l,D,cl 
Cas 



LT 
27rr 



^Res,^Jr(z)(2e)- 



s=0 



LT 



^1;T,0 



4 / 
Z + 1) 



(2z + 7)(2z-3) 



1 



(2z + l)(2z + 3) 



1 



48 
2z-3 



(S+ 1)^+2 

1 4z + 3 



24 



12 (s + 1)^+3 48 (s + 1)^+5 



768V27rr£2 



r^Res.^i 



'2^ + 1^ 



-(2z-3)(2z-l)(2z + 5)Cfl(2 + - 



(5 + 1) 



2+4 



LT 



re 2 



-Ck(3) 



LT 



re 2 



-Cfl(3). 



(41) 



For the term Ai, we have taken the residue at z = 1/2 before taking the sum with respect to s, which give the value 
zero because of the factor (2z — 1). If we sum over s first before taking the residue, then we have the term (fi{z + 1/2) 
which has a pole at z = 1/2. Using the fact that lim (2z — l)Cjf {z + 1/2) = 2, we find that we should add 



LT 



32^/2r£3 



Cas 



In the former case where we do not sum over s before taking the residue, we find that in the high 



to 

temperature region, the Casimir energy has asymptotic behavior 



E, 



D,cl 
Cas 



LT 



and hence the Casimir force has asymptotic behavior 

^D,c\ 3LT 



Cfl(3) (l + | + ...); 



F, 



Cas 



32^2; 



re 2 



C.,(3)(l + ^ 



But we sum over s first, then 



F, 



D,cl 
Cas 



3LT 



re 2 



-Cfl(3) 1 + e -r 



1 



12 6Ci?(3) 



(42) 



(43) 



The term -— = —0.0533 is negative. In 42|, Weber and Gies used worldline numerical method to show that 

12 6Cfl(3) ^ ^ J' 

the first order correction has the same sign as the leading term. Therefore, the result (j42|) is more agreeable. In fact, 

one would also obtain (1421) if one approximate the summation over m by an integral, as in the work of Bordag ^19| 

for the zero temperature Casimir energy. 

For the case of Neumann boundary conditions, using ([29)) . one finds that 



jpl.AT.cl T7.1,-D,cl 
Cas Cas 27^7- 



^Res^ 



r(z)(2£)- 



s=0 



--E. 



l,D,cl 
Cas 



LT 



^Res,^Jr(z)(2£)- 



&2;T.0{z) _ ^1:T.0(2^ + 1) I 1 

r(^) r(22±3) /(, + i)-+^ 

1 f2z+l\ 



<2z+ 



s=0 



(s + 1)" 



l,D,cl 
Cas 



LT 



00 



8V2ret ^ + : 
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which is divergent. If we sum over s first, then 




tC«(3) + — ^ + (l) H2e) + ( i ) - ^^(2) - 2^(1) 



64V2r£2 32\/2re2 8\/2re2 V V2/ 2 V2/ 2 

^^^('21ne+^+81n2-2C+iV 



16V^re2 V 4 2 

where C is the Euler constant. We see that there is a term of order e Ine for the first order correction. For the Casimir 
force, one finds that 



^Caf--7T^CH(3)(l+e(- 



3LT /I 1 



From the discussion above, we see that there is some ambiguities regarding whether one should sum over s first 
before taking the residue. The summation over s appears because we have expanded the logarithm. The divergence in 
s is possibly due to the fact that we have expanded each term in the expansion perturbatively, by treating e and rii /m 
as small variables. When we obtain a divergent sum in s, we need to question whether the perturbation is legitimate. 
For the leading order term of the Casimir force in the medium temperature region, we have found that correct result 
is obtained by summing over s first, performing analytic continuation before taking the residue. However, for the 
first correction term to the Casimir energy in the high temperature region, this seems to lead to a result that does 
not agree with that obtained using numerical approach. The latter might also be special because we are facing the 
divergence series X]^o ^/(^ + which is formally Ci?(l)i infinity. 

V. CONCLUSION 

In this article, we have studied analytically the exact Casimir interaction between a cylinder and a plate at finite 
temperature for a scalar field with Dirichlet and Neumann boundary conditions. By taking the sum of the Dirichlet and 
Neumann Casimir interactions, one obtains the Casimir interaction for electromagnetic field with perfect conductor 
conditions. We rederive the proximity force approximations to the Casimir energy and the Casimir force in three 
temperature regions: the low temperature region where aT <C rT <^ 1, the medium temperature region where 
aT <C 1 <C rT and the high temperature region where 1 <C aT <^ rT. These were compared to the leading terms 
computed from the exact expression. It is shown that in the zero temperature, the medium temperature and the high 
temperature regions, the leading terms derived from the exact expression agree completely with the proximity force 
approximations. In fact, for the zero temperature case, the agreement between the proximity force approximation 
and the exact result has been shown in [l9|. Using a similar perturbation developed in |19l] . we apply the inverse 
Mellin transform representation of the exponential function to derive the first^ analytic corrections to the proximity 
force approximations. In the zero temperature limit, the results obtained in [l9| were recovered. 
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